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ABSTRACT 


Using  randomly  censored  data,  we  develop  a  test  of  the  null  hypothesis 
that  a  new  item  has  stochastically  the  same  residual  lifelength  as  does  a 
used  item  of  specified  age  t^,  versus  the  alternative  hypothesis  that  a  new 
item  has  stochastically  greater  residual  lifelength  than  does  a  used  item 
of  age  t^.  tfe  also  compare  our  test  with  a  related  test,  developed  for  a 
complete-data  model,  in  order  to  study  the  loss  in  efficiency  because  of 
censoring. 


•  •••.*  on  For 


1 

l 


1 


1.  INTRODUCTION 

Suppose  a  cancer  specialist  believes  that  a  patient  newly  diagnosed  as 
having  a  certain  type  of  cancer  has  a  smaller  chance  of  survival  than  does  a 
patient  who  has  survived  5  years  following  a  similar  initial  diagnosis.  How 
can  such  a  claim  be  tested?  To  address  such  a  testing  problem,  Hollander, 
Park,  and  Proschan  (1983)  introduced  the  class  of  "new  better  than  used 
at  tp"  life  distributions  and  the  dual  class  G2  of  "new  worse  than  used  at 
tp"  life  distributions. 

Gj,  the  new  better  than  used  at  tp  alase .  Let  tp>0.  A  life  distribution 
F  (i.e.,  a  distribution  such  that  F(x)  ■  0  for  all  x  s  0)  is  new  better  than 
used  at  tQ  if 

P(x*tQ)  s  F(x)F(tQ)  for  all  x  2  0,  (1.1) 

where  ¥  =  1-F  denotes  the  survival  function. 

G2,  the  new  worse  than  used  at  tQ  olase.  A  life  distribution  F  is  new 
worse  than  used  at  tQ  if  the  first  inequality  in  (1.1)  is  reversed. 

We  now  describe  another  context  in  which  the  concept  of  new  worse  than 
used  at  tp  finds  practical  application.  For  many  electronic  components, 
experience  shows  that  the  failure  rate  of  a  component  is  higher  during  the 
so-called  "infant-mortality  phase"  (say  [0,tp])  than  after  tp.  This  could 
be  the  consequence  of  the  fact  that  there  are  really  two  populations  which 
have  been  mixed  during  manufacture  -  one  population  consists  of  well  made 
units,  while  the  second  population  consists  of  poorly  made  units  whose 
defects  show  up  immediately  upon  initial  use  or  shortly  thereafter. 

The  practical  procedure  used  by  manufacturers  of  such  components  is  to 
"burn-in"  a  portion  of  their  product  during  [0,  tp];  surviving  components  of 
age  tp  sell  at  a  higher  price  than  new  untested  components.  The  purchaser 
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who  needs  higher  than  usual  quality  for  say,  space  vehicle  assembly,  safety 
equipment,  surgical  devices,  etc.,  is  willing  to  pay  for  the  statistically 
higher  quality  units  of  age  tg. 

The  manufacturer  or  one  of  his  high  volume  customers  may  be  highly 
motivated  to  carry  out  the  statistical  test  proposed  in  this  paper  to 
corroborate  his  engineering  experience  and  judgement. 

Gj [G^l  is  related  to  but  contains  and  is  much  larger  than  the  class 
HiPy  of  "new  better  than  used"  ["new  worse  than  used"]  distributions 
defined  below. 

Hj,  the  new  better  than  used  olaee.  A  life  distribution  F  is  new  better 
than  used  if 

F(x+t)  £  P(x)P(t)  for  all  x,  t  i  0.  (1.2) 

H2,  the  new  worse  than  used  olaee .  A  life  distribution  F  is  new  worse 
than  used  if  it  satisfies  (1.2)  with  the  first  inequality  reversed. 

Thus  the  Gj  property  states  that  a  used  item  of  age  tg  has  stochastically 
smaller  residual  lifelength  than  does  a  new  item  whereas  the  Wj  property 
states  that  a  used  item  of  any  age  has  stochastically  smaller  residual  life- 
length  than  does  a  new  item.  Analogous  interpretations  hold  for  G2  and 

The  only  members  in  H^nH2  are  the  exponential  distributions.  Theorem  2 
of  Marsaglia  and  Tubilla  (1975)  shows  that  the  only  members  of  GjnG2  are 
(a)  exponential  distributions,  (b)  survival  functions  F  for  which  F(0)*1 
and  P(tg)=0,  and  (c)  survival  functions  of  the  form  F*(x)=P(x)  for  0sx<tg, 

»  Fj(t0)F(x-jt0)  for  jt0sx<(j«-l)t0,  j=0,l,2,  ....  where  F  is  a  survival 
function  defined  for  x*0.  If  F  has  a  density  function  on  [0,tQ],  then  the 
failure  rate  of  F*  is  periodic  with  period  tg, 
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Some  examples  of  new  better  than  used  life  distributions  are  the  Weibull 
where  F0(t)  «  1-exp (- (It) 0) ,  t2  0,  X2O,  621  and  the  gamma  where 
F0(t)  «  /^x0"1  exp(-Xx)/r(0),  t2  0,  X2  0,  621. 

Let  C*  be  the  class  of  life  distributions  which  are  not  new  better  than 
used  but  are  new  better  than  used  at  t^.  Hollander,  Park,  and  Proschan  (1983) 
give  the  following  method  of  constructing  some  distribution  functions  in  C* . 
Suppose  that  G  is  a  new  better  than  used  distribution  with  failure  rate  function 

rgOO  >  0  for  0  s  x  <  -.  Let  F  have  a  failure  rate  function  rp  satisfying 

(i)rp(x)  s  *G(x)  for  0  £  x  s  tQ,  (ii)rF(x)  »  rG(x)  for  tQ  £  x  <  ®,  and 
(iii)rp(x)  is  strictly  decreasing  for  0  s  t  s  tj,  where  0  <  tj  <  tQ.  Then 
F  c  C*.  To  develop  an  example  of  this  construction,  let  r^(x)  “  1  for 
0  £  x  <  ®,  and  let  rp(x)  *  l-(e/tQ)x  for  0  £  x  <  tQ  and  0  <  6  £  1,  and 
Xp(x)  ■  1  for  tfl  S  x  <  «.  (We  do  not  let  9  exceed  1  since  we  want  to  ensure 
that  rp(x)  remains  positive  as  x  tQ.)  Then  rp  satisfies  (i),(ii),  and 
(iii)  and  thus  F  is  in  C* ,  The  corresponding  survival  function  can  be  written 
as 

f0(t)  =  exp[-(t-8(2t0)'1t2)],  0  <  8  s  1,  0  s  t  <  tQ, 

*  exp[-(t-e(2)_1t0}3,  0  <  9  £  1,  t  2  tQ. 

Hollander,  Park,  and  Proschan  (1983)  derived  a  test  of 

HQ:  F  £  Gj  n  62  (1.3) 

versus 

Ha:  F  c  Gj*  (1.4) 

where 

Gj*  =  (F:  F(x+t0)  £  F(x)F(tQ)  for  all 

x  2  0  and  inequality  holds  for  some  x  2  0). 
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The  null  hypothesis  Hq  asserts  that  a  new  item  has  stochastically  the  same 
residual  lifelength  as  does  a  used  item  of  age  t^.  (Equivalently  F  satisfies 
Hq  if  it  satisfies  (1.1)  where  the  first  inequality  is  replaced  by  an  equality.) 
The  alternative  hypothesis  asserts  that  a  new  item  has  stochastically 
greater  residual  lifelength  than  does  a  used  item  of  age  tQ.  The  test  proposed 
by  Hollander,  Park,  and  Proschan  was  for  the  model  where  we  obtain  a  complete 
random  sample  from  the  life  distribution  F.  In  the  present  paper  we  propose 
a  test  of  Hq  versus  H^  based  on  a  randomly  right- censored  sample.  The  test  is 
derived  in  Section  2.  In  Section  3  we  compare  the  tests  for  the  uncensored 
and  censored  models  and  obtain  a  measure  of  the  loss  in  efficiency  incurred 
because  of  censoring.  Section  4  contains  an  example. 

2. A  TEST  OF  HQ  VERSUS  HA  USING  INCOMPLETE  DATA  ^ 

Let  Xj,  ....  Xn  be  independent  and  identically  distributed  according  to 
a  continuous  life  distribution  F,  and  Yj,  ...,  Yn  be  independent  and  identi¬ 
cally  distributed  according  to  a  continuous  censoring  distribution  H.  Also, 
let  the  X's  be  independent  of  the  Y's.  The  censoring  distribution  H  is 
typically,  though  not  necessarily,  unknown  and  is  treated  as  a  nuisance 
parameter.  In  the  randomly  censored  model  we  do  not  observe  X.,  ...,  X^, 
but  instead,  we  observe  the  pairs  (Z^  6j) ,  i  ®  1,  ...,  n,  where 

Zi  =»  min(X.,Yi) 

1  if  Z.  «  X. 

i  x 

0  if  Z.  »  Y. 

x  i  . 


and 
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Our  test  is  based  on  an  estimator  of  the  parameter 

CD 

T(F)  »  /<F(x*t0)  -  f(x)F(t0))dF(x) 


*  /F(x*t0)dF(x)  -  7FCV  • 


Under  Hg,  T(F)  ■  0  whereas  under  HA,  T(F)  £  0.  In  fact,  since  F  is  assumed 
continuous,  under  H^,  T(F)  is  strictly  less  than  0.  T(F)  gives  a  measure  of 
the  deviation  of  F  from  HQ  towards  H^.  It  is  natural  to  base  a  test  of  HQ 

A  A 

versus  H  on  a  consistent  estimator  of  T(F),  and  we  utilize  T(F  )  where  F 

t\  II  It 

is  the  Kaplan-Meier  (19S8)  estimator.  (Such  an  approach  has  been  used  by  many 
authors  including  Efron  (1967)  in  the  context  of  a  two-sample  location  test, 
by  Koziol  and  Green  (1976)  and  Csorgo  and  Horvath  (1981)  in  testing  goodness- 
of-fit,  and  by  Chen,  Hollander,  and  Langberg  (1983a)  in  testing  whether  new 
is  better  than  used.) 

Under  the  assumption  that  F  is  continuous,  the  Kaplan-Meier  estimator  can 
be  expressed  as 

«  6 
F  (t)  -  S  { (n-i)  (n-i»l)-1)  (i),  t  €  [0,  ZfJ, 
n  U:Z(i)st}  W 


where  =  0  <  Z,,%  <  ...  <  Z,  .  denote  the  ordered  Z's  and  6,..  is  the  6 
(0)  (1)  (n)  (i) 

corresponding  to  Z^j.  Here  we  treat  Z^  as  an  uncensored  observation, 
whether  it  is  uncensored  or  censored.  When  censored  observations  are  tied 


with  uncensored  observations,  our  convention  for  ordering  the  Z's  is  to  treat 
the  uncensored  observations  as  preceding  the  censored  observations. 

Weak  convergence  of  the  Kaplan-Meier  estimator,  regarded  as  a  stochastic 
process,  has  been  established  by  Efron  (1967),  Breslow  and  Crowley  (1974), 
Meier  (1975),  and  Gill  (1983).  Strong  consistency  of  the  estimator  was 
proved  by  Peterson  (1977)  and  Langberg,  Proschan,  and  Quinzi  (1981).  Exact 
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n 

i 
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small -sample  moments  of  Pn,  under  a  model  of  proportional  hazards,  were 
obtained  by  Chen,  Hollander,  and  Langberg  (1982). 

Our  test  statistic  is 
_  def. 

Tn  '  Ttg  = 

For  computational  purposes  we  may  rewrite  TC  as 

n 

T*«  £  n  {(n-k)(n-k*l)_1)4(,t) 

i-1  ||_<k:Z(k)sZ(i)+t0>  J 

iM(n-r) (n-r+l)”1)^(r)j|  •  |^l-{(n-i) (n-m)”1}^! 

£ 

-4  n  {(n-kHn-k+l)"1)  (k). 

{k:Z(k)£V 

i  1/2  c 

Asymptotic  normality  of  n  (T _-T(c)}  can  be  established  under  the 


assumptions 


(A. 1)  The  support  of  both  F  and  H  is  [0,®) . 

(A. 2)  sup{[F(x)]1'E[H(x)]“1,  x  e  [0,®)}  <  ®  for  some  e  2  0. 


Condition  (A. 2)  restricts  the  amount  of  censoring.  For  example,  in  the 
proportional  hazards  model  where  H(t)  »  (F(t) }®,  (A. 2)  implies  that  B  <  1, 
which  means  that  the  expected  proportion  of  censored  observations  $/B*l 
must  be  less  than  .5. 

We  now  state  the  main  result  of  this  section  and  then  sketch  its  proof. 

THEOREM  1.  Assume  that  (A.l)  and  (A. 2)  hold.  Then  n^2{TE  -  T(F)}  converges 

2 

in  distribution  to  a  normal  random  variable  with  mean  0  and  variance  a ■  given 
by  equation  (2.1)  below. 


Sketch  of  Proof.  For  n  ■  1,  2,  . ..,  the  expression  for  -  T(F)  can  be 
written  as: 

T*  -  T(F)  -  /{Pn(x+t0)-F(x+t0)}dFn(x)-/{Fn(x*t0)-F(x+t0)}dF(x) 

♦  / F(x^t0)dFn(x)  -  /F(x*t0)dF(x) 

♦  JPn(x*t0)dF(x)  -  /F(x+t0)dF(x) 

'  ¥n<*0>  *  3*  <V  • 

Unless  otherwise  specified,  all  integrals  range  over  (0,  ®) .  Upon  integration 
by  parts  and  change  of  variables,  we  have: 

'  /F(x*t0)dFn(x)-/F(xn0)dF(x)  =  -/ln(x-t0)dF(x)+/F(x-t0)dFCx),  n  =  1,  2 . 

Thus  for  n  =  1,  2,  .... 

»1/2CT'-T(F))  •  Bnil  ♦  Bn>2. 

where 

Bnl  -  Jn1/2{Pn(x*t0)-F(x+t0)}dFn(x)-/n1/2{Fn(x>t0)-F(x+t0))dF(x), 

and 

“n.2  ’  Inl/2lP„C«*V-f<x*V-lln(,c-to)-p(x-V)-J(f„(V-p(to,,ldFW' 

Using  Theorem  2.1  of  Gill  (1983)  and  arguments  similar  to  those  in  the 
proofs  of  Lemmas  2.2  and  2.3  of  Chen,  Hollander,  and  Langberg  (1983a),  we 
can  show  that 

(a)  B  .  converges  in  probability  to  zero, 

n,  i 

(b)  B  ,  converges  in  distribution  to  the  random  variable 

n,  i 

ZQ  •  t0)-4>(x-t0)  -^(t0))dF(x),  and 


(c)  The  random  variable  ZQ  is  normal  with  mean  0  and  finite  variance 

o2c  =  J/E{<Kx+t0)-<Kx-t0)  -  ^(t0)>  *  (♦(u*t0)-*(u-t0)  -  -^(t0)}dF(x)dF(u).  (2.1) 

Here  (t(t),  te(0,®)}  is  the  Gaussian  process  which  is  the  limit  of  the 
Kaplan-Meier  estimator  regarded  as  a  stochastic  process.  The  mean  of  (t) 
is  zero  and  its  covariance  kernel  is 


E*(t)<Hs)  -  F(t)F(s)  /  (K(z)F(z)  }_1dF(z) ,  0S5<t<«, 

0 


(2.2) 


where  K(t)  =  F(t)H(t).  This  concludes  the  sketch. 

1/2  c 

The  null  asynptotic  mean  of  n  T^  is  zero,  independent  of  the  distri- 

2  £ 

butions  F  and  H.  However,  the  null  asymptotic  variance  of  n  Tn  depends 
on  both  F  and  H  and  thus  it  must  be  estimated  from  the  incomplete  observa¬ 
tions  (Zj.fij),  ....  (Zn,6n).  Under  HQ  it  can  be  shown,  after  straight- 

2 

forward  but  tedious  integration  using  the  expressions  for  o  and  E<t>(t)<J>(s) 

1/2  c 

given  in  (2.1)  and  (2.2),  that  the  null  asymptotic  variance  of  n  Tn  is: 


°c0  =  d/4)p  (t0)/F 5(z)(K(z*t0)> dF(z) 

( 

00  00 

♦  (l/4)F2(tn?fF3(z){K(z))’1dF(z)  -  (l/2)F4(t_) jF3(z) (K(z+t_) }_1dF(z) . 

U  0  0 

-  -  2 

If  there  is  no  censoring,  that  is,  if  K(z)  *  F(z)  for  z  e  [0,  <*>) ,  then  ocQ 
reduces  to  (1/12)P (t0)+(l/12)F2(t0)-(l/6)F3(t0) .  This  expression  agrees 
with  the  null  asymptotic  variance  oQ  of  the  statistic  n  T(Gr)  (where  Gr 
is  the  empirical  distribution  function  of  a  random  sample  from  F)  advanced 
by  Hollander,  Park,  and  Proschan  (1983)  for  testing  HQ  in  the  complete  data 
case.  Expression  (2.3)  can  be  simplified,  by  a  change  of  variable  in  the 
first  and  third  terms,  to 


(2.3) 
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00 

°C0  B  a/4){F(t0)>"2/  F3(u) (K(u) >_idF(u) 

t0 


i-l 


(2.4) 


♦  (l/4)F2(tn)/F3(u){K(u)}'1dF(u)  -  (1/2)/  F^u) (K(u) }_1dF(u) . 

0  t 

*0 


,-l 


*3, 


,-l 


To  obtain  our  estimator  of  o^q  we  introduce  some  notation.  Let 

5  ...  5  denote  the  ordered  Z-values.  Let  Kr  denote  the  empirical 

distribution  function  of  the  Z-values.  Thus  nKn(t)  =  "number  of  Z-values  £  t ." 

Since  Tr  =»  0  when  Z(n)  S  *0-  we  will  assume  our  sample  is  such  that  Z^  >  tQ. 

Replacing  F  by  F  ,  K  by  K  .  and  ®  by  Z,  .  in  (2.4)  yields  the  estimate  o 

n  n  (n)  cn 

defined  by  (2.5). 

o\n  -  [(l/4)(Pn(t0)r2  -  (1/2)] 


rJ  j 


,-l, 


{i:VW(i):SZ(n) 


♦  d/4)P2(t  )  J 

u  (i:lsiST (n) } 


l  >  <Rn<W(i)  ->  >"i{Fn(W(i)  )  -  VW(i-l) )  > 


(2.5) 


where  »(0)  =  0  <  »(1)  <  W(2)^< 


<  W,  ,  are  the  ordered  observed 
(T(n)) 


failure  times,  and  -rfn)  *  £  6.  is  the  total  number  of  failures  among  the 

i=l  1 

n  observations. 


•  2 
We  are  unable  to  prove  that  o cn  consistently  estimates  under  our 

a  2 

(A.l),  (A. 2)  assumptions,  but  we  have  investigated  properties  of  ocn  via 

a  2 

a  limited  Monte  Carlo  study.  Table  1  investigates  the  accuracy  of  ocn  as 
an  estimator  of  and  the  accuracy  of  the  normal  approximation  in  the 
cases  where  F  is  exponential  with  scale  parameter  1  and  H  is  exponential 


10  - 


with  scale  parameter  X,  for  the  choices  X».l  and  X=l/3,  with  t0®.6,l,  and 

a 2 

n=100,  150,200.  Column  2  of  Table  1  gives  the  average  value  of  ocn>  averaged 
over  1,000  Monte  Carlo  replications.  Column  3  gives  the  sample  standard 

A  n 

deviation  s  of  the  1,000  values.  It  is  seen  that  o  tends  to  be  below 

cn  cn 

2 

the  true  value  oc^,  but  the  estimator  improves  as  n  increases. 

The  approximate  o-level  test  of  Hq  versus  H^,  which  rejects  Hq  in  favor 

of  if  n1/2T^oc^  s  and  accepts  Hq  otherwise,  is  called  the  NBU-t^  test. 

The  approximate  a- level  test  of  HQ  versus  the  alternative  that  a  new  item  has 

stochastically  smaller  residual  lifelength  than  does  a  used  item  of  age  t^ 

is  called  the  NWU-trt  test.  The  NWU-tn  test  rejects  H_  if  v}^Tco  *  a  z 
0  0  J  0  n  cn  o 

and  accepts  HQ  otherwise.  Here  z^  is  the  upper  a-percentile  point  of  a 

standard  normal  distribution.  Columns  4, 5, 6, 7, 8, 9  pertain  to  the  convergence 

to  asymptotic  normality  of  the  standardized  test  statistic  T*  *  v}^2Tco~l . 

n  cn 

Columns  4,5,6  give  estimated  probabilities  of  the  events  {T*  £  -zo),  and 
columns  7,8,9  give  estimated  probabilities  of  the  events  (T*  a  zo>, 
a  *  .10,. 05,. 01.  It  is  seen  that  the  convergence  to  asymptotic  normality 
is  slow.  The  probability  o  assigned  to  the  event  {T*  £  -zq}  by  the  normal 

A 

approximation  is  less  than  the  corresponding  Monte  Carlo  estimate  P{T*  £  -z^} . 

Thus  the  NBU-tQ  test  based  on  the  normal  approximation  tends  to  give  P  values 

that  are  less  than  the  true  P  values.  The  probability  a  assigned  to  the 

event  {T*  a  zq}  by  the  normal  approximation  is  greater  than  the  corresponding 

Monte  Carlo  estimate  P{T*  a  z  }.  Thus  the  NWU-trt  test  based  on  the  normal 

a  u 

approximation  tends  to  give  P  values  that  are  greater  than  the  true  P  values. 
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3.  EFFICIENCY  LOSS  DUE  TO  CENSORING 


In  this  section  we  study  the  efficiency  loss  from  censoring  by  comparing 

A 

the  efficacy  of  the  test  based  on  Tr  =  T(Gn)  for  the  uncensored  model  with 

g 

the  efficacy  of  the  test  based  on  Tr  for  the  randomly  censored  model.  Also 
we  present  a  monotcnicity  property  of  the  efficiency  as  the  amount  of  censor¬ 
ing  increases. 

Since  Tn  and  Tn  have  the  same  asymptotic  means,  the  value  of  1-R,  where 

R  =  ‘f.X-  V  ■ 


can  be  taken  as  a  measure  of  the  efficiency  loss  due  to  censoring  (cf.  Chen, 

2  2 

Hollander,  and  Langberg,  1983a).  Here  oQ  and  ocQ  are  the  null  asymptotic 
variances  of  n*^2Tn  and  respectively. 

We  consider  the  case  where  the  censoring  distribution  is  exponential 
with  parameter  X  and  the  life  distribution  is  exponential  with  parameter  1. 
To  satisfy  Condition  (A. 2)  we  take  X  <  1.  Then  we  have 


?  -tn  -2tn  -3t 

at.  «  (l/12)e  u  ♦  (l/12)e  u  -  (l/6)e  u, 


°c0  *  U/4(3“x»  *  (e 


-(l-X)t. 


■2tr 


+  e 


-  2e 


"(3-X)t, 


and 


R  *  (l-(X/3) }(l+e  -2e 


_2to  Ato  -to 

°)/(e  °*e  -2e  1 


(3.1) 


In  Table  2,  values  of  R  given  by  (3.1)  are  presented  for  several 
choices  of  t^  and  X  <  1.  The  table  shows  that  as  X  decreases,  the  value  of 
R  increases  to  1.  Note  that  X  ■  0  implies  no  censoring.  The  table  also 
shows  that  as  tQ  increases,  the  efficiency  of  T^  with  respect  to  Tfi 


decreases. 


The  next  theorem  shows  that  if  is  stochastically  smaller  than 
(i.e.,  there  tend  to  be  more  censored  observations  with  censoring  distri¬ 
bution  than  with  Hj),  then  the  efficiency  of  with  respect  to  Tn  under 
Hj  is  greater  than  under  H2> 

THEOREM  2.  Assume  that  Hj  £  H2,  where  Hj  and  H2  are  censoring  distributions 
Then  Rj  a  R2,  where  Rj  =  ep  ^(Tn,  Tn5  and  R2  “  eF,H2(Tn’  V* 

Proof.  Since  the  numerator  of  e_  u  (Tc,  T  )  i  *  1,2,  does  not  depend  on 

r,rK  n  n 

H,  it  suffices  to  show  that  S  °cO^H2^  wbere  °c0^  *S  g*ven  by  C2*4) 

(Recall  that  K  *  FH.)  Since  H ^  £  H2  we  have  [Hj]"1  £  Note  that 

(H)  is  increasing  in  [H]"1.  Thus,  it  follows  that  ^o^P  5  °cO^H2^  * 
Consequently,  the  desired  result  follows. 

From  Theorem  2,  it  immediately  follows  that  the  maximum  value  of  R  is 
equal  to  1  and  is  achieved  when  H(x)  2  1  for  x  a  0,  that  is,  when  there  is 
no  censoring. 


Table  2.  Efficiency  of  T^  with  respect  to  Tr  when  F(x)  *  e”x,  x  i  0 

and  H(x)  -  e”*x,  x  £  0,  X  >  0 


1/2 

1/4 

1 

1/8 

1/20 

1/60 

1/100 

.9008 

.9514 

.9759 

.9904 

.9968 

.9981 

.7265 

.8583 

.9280 

.9709 

.9903 

.9942 

.5823 

.7711 

.8804 

.9509 

.9835 

.9901 

.4652 

.6907 

.8338 

.9307 

.9765 

.9858 

.3325 

.5843 

.7673 

.9003 

.9657 

.9793 

.1931 

.4442 

.6686 

.8520 

.9482 

.9686 

.1146 

.3411 

.5854 

.8077 

.9314 

.9583 

.0688 

.2639 

.5146 

.7670 

.9154 

.9484 

4.  AN  EXAMPLE. 


The  data  in  Table  3  are  an  updated  version  of  data  analyzed  by  Koziol 
and  Green  (1976) .  The  data  correspond  to  211  State  IV  prostate  cancer 
patients  treated  with  estrogen  in  a  study  by  the  Veterans  Administration 
Cooperative  Urological  Research  Group  (1967).  By  the  March  1977  closing 
date,  90  patients  had  died  of  prostate  cancer,  105  had  died  of  other  diseases 
and  16  were  still  alive.  The  latter  121  observations  will  be  treated  as 
censored  observations  (withdrawals) . 


Table  3.  Survival  times  and  withdrawal  times  in  months  for  211  patients 
(with  number  of  ties  given  in  parentheses) 


Survival  times:  0(3),  2,  3,  4,  6,  7(2),  8,  9(2),  11(3),  12(3),  15(2),  16(3), 
17(2),  18,  19(2),  20,  21,  22(2),  23,  24,  25(2),  26(3),  27(2),  28(2),  29(2), 
30,  31,  32(3).  33(2),  34,  35,  36.  37(2),  38,  40,  41(2),  42(2),  43,  45(3),  46, 
47(2),  48(2),  51,  53(2),  54(2),  57,  60,  61,  62(2),  67,  69,  87,  97(2),  100, 
145,  158. 


Withdrawal  times:  0(6),  1(5),  2(4),  3(3),  4,  6(5), 7(5),  8,  9(2),  10,  11, 

|  12(3),  13(3),  14(2),  15(2),  16,  17(2),  18(2),  19(3),  21,  23,  25,  27,  28,  31, 

32,  34,  35,  37,  38(4),  39(2),  44(3),  46,  47,  48,  49,  50,  53(2),  55,  56,  59, 
61,  62,  65,  66(2),  72(2),  74,  78,  79,  81,  89,  93,  99,  102,  104(2),  106,  109, 
119(2),  125,  127,  129,  131,  133(2),  135,  136(2),  138,  141,  142,  143,  144, 
148,  160,  164(3). 


Koziol  and  Green  (1976)  stated  that  experience  had  suggested  that  had 
the  patients  not  been  treated  with  estrogen,  their  survival  distribution  F 
for  deaths  from  cancer  of  the  prostate  could  be  taken  to  be  exponential  with 
a  mean  of  100  months.  With  this  in  mind,  various  authors  have  applied 
goodness-of-fit  tests  to  the  prostate  cancer  data.  A  recent  reference  is 
Csorgo  and  Horvath  (1981)  whose  procedures  indicate  lack  of  support  for  the 
simple  null  hypothesis  that  F  is  exponential  with  mean  100.  (Csorgo  and 
Horvath  also  discuss  other  references  relating  to  tests  of  this  simple  null 
hypothesis.) 


Tests  of  the  composite  null  hypothesis  of  exponent iality  (with  unspeci¬ 
fied  mean)  were  performed  by  Chen  (1981),  Chen,  Hollander,  and  Langberg 
(1983a),  and  Chen,  Hollander,  and  Langberg  (1983b).  These  tests  indicated 
lack  of  support  for  the  hypothesis  of  exponential  aging.  Chen,  Hollander, 
and  Langberg  (1983b)  also  plotted  an  empirical  mean  residual  life  function 
for  the  data  of  Table  3.  Their  plot  tends  to  decrease  up  to  around  25  months, 
then  tends  to  increase  up  to  about  70  months,  and  then  decreases  again. 

The  null  hypothesis  Hq  (1.3)  of  this  paper  is  appropriate  if  one  has 
a  priori  reasons  to  suspect  that  a  patient  after  tQ  months  would  have  sto¬ 
chastically  the  same  residual  lifelength  as  a  new  patient.  However,  we  are 
not  aware  of  any  such  a  priori  notions  in  this  prostate  cancer  setting,  and 
there  is  no  natural  value  of  tQ.  Hence,  our  test  is  performed  primarily  for 
purposes  of  illustration.  With  the  choice  tQ  »  60,  we  find  T^jj  *  .0116, 
o*211  *»  .3354,  and  (211)  1/,2t2h°c211  *  ,290»  8  value  supporting  HQ. 
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